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Abstract 

Spin transport driven by an external electric field in uniform metallic ferromagnets with the 
spin-orbit interaction arising from random impurities is studied microscopically. Spin relaxation 
torque T is shown to be written by spatial derivatives of the electric field, but with anisotropy 
arising from the magnetization. The field-driven contribution of the spin current is also anisotropic. 
The diffusive spin current is shown to be written as a gradient of the spin chemical potential, and 
the linear-response expression for the spin chemical potential is derived. It is discussed that the (5 
term in the spin transfer torque can also be anisotropic. 

PACS numbers: 72.10.Bg 72.25.Rb 
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I. INTRODUCTION 



Spintronics aims at using the information carried by the electrons' spin in sohds. For 
this purpose, estabhshing rehable methods to create, transfer and detect the spin current 
is an urgent task. Compared to the charge transports, spin transports have one serious 
fundamental difficulty. That is the non-conservation of the spin in solids. This limits the 
range of the spin transmission to be less than the spin diffusion length, which is typically 
yum scale in metals. 

The non-conservation of the spins is expressed by a source term in the continuity equation 
for the spin 

s" + v-jr = r". (1) 

Here s and js are the spin density and spin current, respectively, a = x,y, z is the spin 
direction and T is the spin relaxation torque resulting in the non-conservation of the spin. 
In the most cases in metals, the dominant origin of T is the spin-orbit interaction. 

Although the relaxation torque term is essential in spin transports, it has so far been 
treated only on the phenomenological ground. The continuity equation is equivalent to 
the Boltzmann equation, which is useful in discussing the spin transports. The Boltzmann 
equation for the distribution function of each spin channel was discussed by Son et al. [l| 
and later by Valet and Fert ^ in the context of the giant magnetoresistance in multilayer 
systems. In their analysis, they approximated the spin relaxation torque as proportional to 
the inverse of a spin relaxation time Tgf and to some unknown function representing a driving 
force for the spin accumulation. The driving force was written in terms of what they called 
the spin chemical potential fig- The relaxation torque was approximated as = fig/Tsi- 
They argued that /is satisfies the diffusion equation, V^/is = — £^^yUs, with the diffusion 
length 4f oc ^/^s^■ Microscopic calculation for /ig has not been done so far. 

The diffusion equation for the spin has been widely used to discuss recent spin transports 
in metallic junctions jsl. The decay of spin transport has been confirmed in non-local 



spin injection experiments 



-6|, which indicate that the spin diffusion decays with a decay 
length of 350-500nm in Cu and lOOnm in Au at room temperature. Although the spin 
diffusion equation appears to be so far successful, the phenomenological treatment of the 
spin relaxation term and the spin chemical potential must be improved to consider the spin 
transport seriously. 



Besides diffusive spin current, there is another spin current that is driven by an effective 
field. In contrast to the diffusive one, this field-driven contribution should not decay in 
uniform (single domain) ferromagnets, since the ratio of the spin current and the charge 
current is determined by the spin polarization ratio of the material, which is a statistical 
mechanical quantity. The field-driven (local) spin current and the diffusive spin current 
behave differently, as was recently demonstrated theoretically in the case of the inverse spin 
Hall effect Q- 

In the field of the current-driven magnetization dynamics, the spin relaxation torque has 



been studied from the microscopic viewpoint jsHlOf. In this context, Eq. ^ gives the 



expression for the torque acting on the spin density s as 

= -V • + T". (2) 

In the adiabatic limit, i.e., slowly varying magnetization, and under uniform current, the 
:irst term reduces to V ■ = (P/2e)(j ■ V)s°, where P is the spin polarization of the current 
9], [lo|, namely to the adiabatic spin-transfer torque. When the spin-relaxation sets in, the 



conduction electron no longer follows the magnetization profile, and new contribution to the 
torque arises from the T term. This torque was shown to be 

T=-P^{sx{j.V)s), (3) 

where /3 is a coefficient inversely proportional to the spin relaxation time Tg js], [lo|. This 
torque, called /3 term, turned out to be essential in determining the efficiency of the current- 



driven domain wall motion 
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13| . The magnitude of the parameter /3 has recently been 
intensively studied experimentally by measuring the domain wall speed under current [l^ 



151. Theoretical formulation for estimating B in the first-principles calculations was carried 

n 

out recently |16l |. 

The spin relaxation torque has been studied also from the viewpoint of how to define the 
spin current. It was discussed that the spin relaxation torque contains a term written as a 

n 

divergence of the torque dipole density, P p/7| • Generalized argument was given by Shi et 



al. 



18|, where they discussed that the z component of the relaxation torque is written as a 



divergence, 

= -V ■ P, (4) 
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if the system has the inversion symmetry. This means that the total torque integrated over 
the system should vanish. Shi et al. also argued that if the relaxation torque is a divergence 
of P, one can define a spin current that is conserved. In fact, defining js = js + P, the 
continuity equation ([T]) reduces to + V ■ js = 0. The explicit form of the torque dipole 



density was not calculated in Ref. [18|. Obviously, in the presence of the inhomogeneity of 
the magnetization, the /3 torque (Eq. Q) cannot be written as a divergence, and thus it 
indeed represents the spin angular momentum lost by the spin relaxation. 

The result that the spin relaxation torque is given by a derivative of the applied electric 
field E is understood as follows. The spin relaxation torque should of course vanish when 
E = 0. It cannot be directly proportional to E, since field-driven spin current in uniform 
ferromagnets should not decay. Therefore, the simplest expression for the relaxation torque 
is a derivative of the field. It is not, however, obvious whether it should be always written 
in a rotationally invariant way, or if it can be anisotropic, since the rotational invariance is 
broken in uniform ferromagnets because of the magnetization. 

The first aim of the present paper is to calculate the spin relaxation torque microscopically 
in the presence of the applied electric field. The spin relaxation mechanism we take into 
account is the spin-orbit interaction due to random impurities. Our explicit calculation 
reveals that the spin relaxation torque is not always rotationally symmetric in uniform 
ferromagnets, but is generally given by 

r' = ^i\/-E) + S^id,E,), (5) 

where z axis is along the magnetization, 7 and ^7 are coefficient proportional to the inverse 
spin relaxation time. The spin relaxation torque is, therefore, anisotropic. Relaxation torque 
of Eq. ([5]) indicates that the torque dipole density is given by 

P = -{-fE + 6-f{n-E)n), (6) 

where n represents the direction of the magnetization. These anisotropic behaviors of the 
transport quantities is common when spin-orbit interaction exists, as is well-known in charge 
transport as the anisotropic magnetoresistance (AMR) [l^. 

The second aim of the paper is to study the spin current on the same microscopic footing 
as the relaxation torque. We show that the spin current is made up of a field-driven contri- 
bution which is local and the diffusive one with nonlocality. The field-driven contribution 



is anisotropic like the AMR effect for the charge current (we call the effect as spin AMR 
effect). The diffusive contribution is given as a gradient of a spin chemical potential, /ig. We 
will derive the linear-response expression for the spin chemical potential. Our microscopic 
study on the spin current demonstrates the validity of the half-phenomenological treatments 

a. 

II. MODEL 

We consider the conduction electron system taking account of the spin-orbit interaction, 
the impurity scattering without spin flip, and the applied electric field. The Hamiltonian of 
the system is given as if = ifo + -f^so + -f^cm + -f^imp; where Hq is the free electron Hamiltonian 
including the uniform magnetization, Hso is the spin-orbit interaction. Hem is the interaction 
with the gauge field representing the applied electric field, and i^imp is the spin-independent 
impurity scattering. The free part reads 

^0 = "^ ekaci^Cka, (7) 

ka 

where the electron creation and annihilation operators are denoted by and c, respectively, 
eka = ^''^F — o'M, ep is the Fermi energy, M is the spin splitting due to the magnetization 
and cr = ± represents the spin. The spin-orbit interaction is represented by -f^so = -^so° + 
Hso"^, where {ak [k = x,y,z) is the Pauli matrix) 

^so° = E '^^f' J dH^^vi'J){c^ V, akc) (8) 

ijk 

Hso^ = -eJ2 ^^ik j d^x{Viv'£>)A,{x, t)(cV,.c). (9) 

ijk 

(We suppress the spin index when obvious, namely, c = (c+,c_). ) The spin-orbit potential 
is assumed to arise from random impurities and to depend on the spin direction [k). 
The averaging over the spin-orbit potential is carried out as 

(^(„^)(p)^W(-p')>i = n,o\j5^,p'k,, (10) 

where tt-so and Ago are the concentration of the spin-orbit impurities and the strength of the 
interaction, respectively. The average of the spin-orbit potential at the linear order is zero in 
our model, and thus we do not take account of the the anomalous Hall and spin Hall effects. 
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We consider a case where electric field E{x,t){= —A{x,t)) is position and time dependent. 
The electromagnetic interaction is written as 



Hem = -— ^iMQ^ ^^)(c^fc^aCfc+|), (11) 

k.q i 



where Q is the frequency of the electric field. We will consider the limit of small Q and small 
q. The scattering by the normal impurities is represented by 

N- 

^ * imp 

^-p = EE|^'^'"''^''''^'^'^^' (12) 

i=l kk' 

where Vi represents the strength of the impurity potential, Ri represents the position of 
random impurities, A^imp is the number of impurities, and = V/a^ is number of sites. 
To estimate physical quantities, we take the random average over impurity positions in a 
standard manner j^. 

To derive the spin continuity equation, Eq. ([T]), we derive the equation of motion for the 
spin density, s{x,t) = (^c\x,t)(Tc{x,t)) (( ) represents the quantum average). The time 
development of the spin density reads 

s"" = i{[H,c^]a''c + c^a''[H,c]) , (13) 

where H is the total Hamiltonian of the system. The commutators are calculated as in 
Appendix lAl and Eq. ( fT3l) turns out to be Eq. ([T]), namely 

= -V-jT + T", 

with the spin current given as 

jr=jf^'"+jr", (14) 

where 

^iS^+ltt, (15) 



and 
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The relaxation torque reads 

= 7;: + 7;^'", (17) 

where 

r,l = tY^ ei,keaik{V,v^^) (c^ai Vj c) , (18) 

ijkl 

^2eJ^ e^,keM'^^vi'J)A, (c^aic) . (19) 

ijkl 

The spin relaxation torque depends on the definition of the spin current. For instance, if 
we redefine the spin current as ji."' = — where G is a vector, the continuity equation 
([1]) becomes s° + V ■ j!," = T'", where the relaxation torque reads T'° = + V ■ G". This 
ambiguity of spin current definition of course does not affect physical quantities such as the 
total torque acting on the spin density, which is given by s. 

III. SPIN RELAXATION TORQUE 

We calculate the spin relaxation torque as a linear response to the applied electric field. 
The uniform magnetization is chosen as along z axis. The spin-orbit interaction is included 
to the second order. The contributions to the relaxation torque, Eq. f[T7|) . are shown in Fig. 



Vso 




FIG. 1: Feynman diagrams representing the relaxation torque. Solid lines represent the electron 
Green's function with the lifetime (to-) included, Vgo represents the spin-orbit interaction (double 
dashed line) and dotted line represents the interaction with the gauge field (A). The first two 
diagrams are the contributions to 7^" and the last diagram is the contribution to %o'°^- The vertex 
marked by cross represents the relaxation torque. 

[U The leading contribution for small l/lepr) and qi {i is the electron mean free path) turns 



7 



out to be the first diagram in Fig. [H which reads (see Appendix [B] for details) 



2_ , 2±_ 

ovr 

J2 k\kr Yl <y9l,.9l>,-.{9l',-.f + C.C., (20) 



457r 



X 

fcfc' o-=± 



where gl.^ and (7^^. are the retarded and advanced electron Green's functions, respectively, 
carrying the wave vector k and spin a with zero frequency. As we see, only z component of 
the torque is finite. The Green's functions include the lifetime arising from the self-energy 
process due to normal impurities and the spin-orbit interaction. The inverse lifetime for the 
electron with spin cr(= ±) is given as 

Ta~^ = 27mivfu^{l + K^^^ + K_l7<7), (21) 

where u^r is the spin-resolved electron density of states, rii and vi are the concentration and 
the potential strength of the impurities, n^^a = kp^ and k± = kp_^_kp_ are 

dimensionless ratios of the spin-orbit interaction to the normal impurity scattering {kp^^ 
is the spin-dependent Fermi wavelength), and 70- = The total relaxation torque is 

therefore given by (Eq. (E])), 

r' = ^{V-E) + 6^{d,E,), (22) 

where 

7 = ^^nso\so^iy+i^^kl^kl_ ^ (^kp^T^ (23) 

cr=± 

and 5'j = 7/3. The parameter 7 is proportional to the spin flip rate due to the spin-orbit 
interaction. Our result indicates that the relaxation torque is zero in uniform ferromagnet 
when uniform electric field is applied. Thus the spin current does not decay in this case. In 
fact, the static solution of Eq. ([1]) with = is jg = constant. 

The degree of the asymmetry, ^7/7, is not universal but is model dependent. For instance. 



in the case of junction with weak electron hopping at point-like leads, ^7 vanishes 20|. 



IV. SPIN CURRENT 



We here calculate the spin current within the same formalism. Within the linear response 
theory, the spin current is calculated by estimating the Feynman diagrams shown in Figs. 
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[2] and [3l which correspond to the field-induced contribution and the effect of the diffusive 
electron motion, respectively. 



t'so 




FIG. 2: Feynman diagrams representing the local contribution to the spin current (the first two 
terms in Eq. (j37p ). The first three diagrams correspond to and the last two diagrams represent 
the contribution from the anomalous velocity, Dotted and double dashed lines denote the 
interaction with the applied electric field and the spin-orbit interaction, respectively. The vertex 
marked by cross represents the spin current. 




FIG. 3: Left: The vertex correction contribution to ji , resulting in the diffusive spin current 
(the last term of Eq. (|37p ). Right: T„ui, which is a ladder process of the successive electron 
scattering by the normal impurity and the spin-orbit interaction (represented by thick dotted 
lines) connecting the spin indices a and a' . 



We first estimate the normal part of the spin current, j^^^'^ (Eq. 0151) ). shown in the 
first two diagrams in Fig. [21 The contribution jf^^'^ is defined including the anomalous 
contribution from the electromagnetic gauge field, jfl^- Its dominant contribution in the 
limit of small Q and small q is calculated as (see Appendix [C] for detail) 

. 1 



.{n),z 
3s,i 



—I 



2t[ rn^ 



kg j a 



5^ A;,A;;.^^_2,,^fc+|,^5fc'-2,^'5fc'+2 ^/riitf r,„,(g, fi) 



(24) 



In Eq. (El]), the first term is the contribution shown in the left of Fig. |2l and the second 
term is the contribution from the vertex correction (Fig. [3]). 
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The factor Tcra'iQ, ^) contains all the vertex corrections due to the normal impurities and 
the spin-orbit interaction shown in Fig. |3l The equation of motion for F^ct' is derived in the 
same manner as in Ref. 2l| carried out in the context of quantum correction (the diffusion 



without the spin-orbit interaction was considered in Ref. 



22|). The equation is obtained as 



F,^^ = (1 + + n^F^^) + K_Ln_<xr_^,^ 
r.,-a = ^±{1 + n_,F„^,_,) + (1 + /€.)n^F^,„,, (25) 



where 



k 

^ [1 - {D^q^T^ + K, + K^'y^)] , (26) 

where Dcr = ^3^2 is the diffusion constant. Neglecting quantities of order of (k^,k±)^, 
Eq. (125!) is solved as 

1 + - (1 + 2k,)U_^ 



^'^^ [l-(l + «:,)n+][l-(l + «:,)n_] 
F = II (27) 

[i-(i + «:,)n+][i-(i + «:,)n-]' ^ ^ 

Using Eq. fl26|) . we obtain T^^/ as (assuming the rotational symmetry for the wave vectors 
when averaging over the spin-orbit potential) 

1 



F 



^^■-"^ [D^qH^ + «:^7+] [D^q^r^ + k^iA ' ^^^^ 

By use of Eq. ( 128|) and summing over the wave vectors, the normal spin current, Eq. (!24|l . 
reads 

= ^'sE. - (29) 

where = e ^^(±)D±z/± is the bare spin conductivity divided by e. The first term of Eq. 
( 129|) is the field-driven contribution. The second gradient term is a diffusive contribution 
(vertex corrections), arising from the spin accumulation. The effective potential describing 
the spin accumulation, /is, reads 

.../A',(.-.0(v.i.)M. (30) 
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where x is a correlation function arising from the electron diffusion, given as {V is the system 
volume) 

Here a±. = eD^v^ is the spin-resolved Boltzmann conductivity divided by e, and the corre- 
lation length is given as 



4,. = a/^D^. (32) 

The lifetime of the spin a electron reads r^^a = T~a / {,i^i.la)- Defining /ig = /i+ — we see 
that spin-resolved effective potential satisfies 

+ (4,.)-')/i. = -(Ta{V ■ E). (33) 
In three-dimensions, the correlation function reads 

^(-) = 4^E(±)-±^"'^"'^"- (34) 

The local part of the spin current arises also from the anomalous spin current due to the 
spin-orbit interaction, defined in Eq. f lT6|) . This contribution is calculated by evaluating the 
last two diagrams in Fig. [2] as (see Appendix [C]) 

3^^ = 6a,{l - 5.,.)E,, (35) 

where 

6a, ^ In^oXso'- y2{±){kp±)\i^±fr^. (36) 
9 m ^ 

This spin-orbit correction to the spin conductivity is anisotropic, resulting in a spin version 
of the anisotropic magnetoresistance (AMR) effect, namely, spin AMR effect. 

From Eqs. (12^ ( l28l) ( l35l) . the leading contribution to the spin current for small q and Q 
is obtained as the sum of the local part driven by the electric field and the diffusive part as 

= a,E - 6a,{n ■E)n-Vfi,, (37) 

where a, = + Sa, and n is the unit vector along the magnetization. In terms of the angle 
9 defined by cos 6' = {n ■ E)/E, the magnitude of the field-driven (local) current reads 



j'r' = \/(^s||)^ + ((crsx)^ - (^s||)2) sin^ 9, (38) 
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where as\\ = cr^ and as± = (Ts- When the degree of the anisotropy is small, the spin current 
becomes 

We define the magnitude of the spin AMR as 

Aps ^ Ps|| - Ps± ^ ^^s/q-s 

Ps± ~ psL 1 - 6a J 
where = (osa)'^ {a =||, ±). 



V. SPIN INJECTION 



We have thus derived the explicit expression for the spin chemical potential within the 
linear response theory. Let us apply Eq. ( 130|) to a ferromagnetic-normal metal junction 
with an insulating barrier, used in the nonlocal spin injection experiments jsl, depicted in 
Fig. m^a). When the voltage is applied perpendicular to the interface (we choose the x axis 
in this direction), the electric field is uniform inside the ferromagnet and the normal metal 
except at the interface. Writing the voltage drop at the interface (chosen as at x = 0) by 
VpN, we obtain 

V-E^ 6{x)VF^/d, (41) 

where d is the width of the interaface, which is treated as small enough compared with 
the electron mean free path, resulting in the delta function in V ■ -E. In totally unpolarized 
non-magnetic metals, namely, if 0"+ = cr_ and -D+ = D_, the correlation function in Eq. (13 ip 
always vanishes. As is naively guessed, therefore, spin injection thus requires an effective spin 
polarization close to the interface, induced by the exchange interaction with the ferromagnet. 
This spin polarization is expected to be localized within a short distance of a few lattice 
constants from the interface. Let us approximate the interface polarization by introducing 
spin-dependent diffusion constant and the density of states, D^^ and respectively, at the 
interface. The long-range behavior of the spin correlation function in the non-magnetic side 
is then obtained as 
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where £s in the spin diffusion length in the normal metal ( 4 is a long (~ fim) length scale 
and thus does not depend on the spin). We therefore obtain from Eq. fl30|) the chemical 
potential as 



4:7l\x\ 



(43) 



where qs = {^^aD^ly(j)eVF^)Ap-!<(/d, is the spin accumulation rate at the interface (per unit 
time), and ApN is the area of the junction. This result of fig is consistent with intuitive and 
phenomenological results of the spin injection in the perpendicular structure shown in Fig. 
m^a). In contrast, when the voltage is applied parallel to the ideal interface as shown in Fig. 
m^b), spin injection does not occur since V ■ -E = at the interaface and thus /xP^ = 0. 



Vfn 



<7s 



iF N 



(b) 



FIG. 4: (a) Creation of diffusive spin current (spin injection) by applying the electric voltage 
perpendicular to the F-N interface, (b) When the voltage is applied parallel to an ideal interface, 
no spin current is induced since V • £^ = 0. 



VI. TOTAL TORQUE AND ASYMMETRIC /3 TERM 

The continuity equation ([1]), indicates that the spin polarization (magnetization) changes 
due to the spin relaxation. (Change of the magnetization magnitude is a feature of the 
itinerant magnetism.) By use of Eqs. (l5i) (l33i) (1371) . we see that Eq. ([T]) results in 

s^ = ^rV-E + S^rV.E, + $^(±)77^, (44) 

-I- l^s,±J 

where 7^ = 7 — Sa^ and ^7^ = ^7 + Sa^. General case with uniform magnetization along any 
unit vector n is given by (s = sn) 

i = n ( 7.V ■ + 57.V||E|| + ^(±)^£±^ j , (45) 
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where E\\ = n ■ E and Vy = n ■ V. 

In addition to the change of the magnitude, Eq. (H5l) . there is a torque, which is per- 
pendicular to n. Such torque arises when the magnetization is not homogeneous, and plays 
important roles in current-induced magnetization dynamics. We have carried out the calcu- 
lation of the current- induced torque done in Ref. [lOj] on the same footing as the derivation 
of Eq. ([5]). As a result, we found that the /3 term becomes asymmetric as (see Appendix [D] 
for details of the calculation) 

Tli'^^'" = X (j ■ V)s + 5Ps X (Jii V|i)s]°, (46) 

where j\\ = n j is the current along the local magnetization and 6(3/P = —1/5 in the present 
model. The spin transfer torque due to the spin-orbit interaction is thus different from that 
due to the spin-flip scattering. The expression of the total torque allowing for the spatially 
varying current density and the magnetization is therefore obtained as 

P P 
s = - — (V • j)n - - [/3n X {j ■ V)n + 5/3n x (jy V||)n)] 

+ n ■ E) + 5^r{d\\E\\) + ^(±)(4,±)-V± j • (47) 

This expression clearly demonstrates that the spin relaxation torque requires some inhomo- 
geneity either of the applied current or the spin structure, in addition to the spin-orbit (or 
spin flip) interaction. Totally homogeneous system does not relax. The last term in Eq. (l47|l 
gives useful information for measuring the spin accumulation induced by the spin current. 



VII. CONCLUSION 



We have carried out a microscopic calculation of the spin relaxation torque and the 
spin current induced in disordered ferromagnetic metals by the applied electric field. The 
spin-orbit interaction arising from the random impurities is included as a source of spin 
relaxation, and inhomogeneity of the applied electric field is taken into account. We found 
that the spin relaxation torque in the uniform magnetization case is written as a divergence 
of the electric field plus an anisotropic term. The spin current was shown to be made up 
of field-driven (local) and diffusive (nonlocal) contributions, the latter written as a gradient 
of a spin chemical potential. We have derived a general linear response expression for the 
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spin chemical potential. The spin injection effect was brieffy discussed based on our results. 
When the analysis is applied to the inhomogeneous magnetization case, we argued that the 
(3 torque in the current-induced magnetization dynamics can be anisotropic. 

Before finishing, we emphasize that the expression for the spin current and are mean- 
ingless without specifying the physical observable to be measured. In the inverse spin Hall 
effect, which was originally proposed as 
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jfj_ oc e^vpjf^jy; it has recently been demon- 



strated that the charge current is not directly proportional to the spin current 
for the spin current only does not therefore provide physical information. 



3, 



25|. Solving 



Acknowledgment 

The authors thank E. Saitoh, J. Shibata, H. Kohno, S. Murakami for valuable discus- 
sions. This work was supported by a Grant-in-Aid for Scientific Research in Priority Areas, 
"Creation and control of spin current" (1948027), the Kurata Memorial Hitachi Science and 
Technology Foundation and the Sumitomo Foundation. 



Appendix A: Derivation of the spin continuity equation 

We calculate the commutators which appears in Eq. ( fT3l) with the total Hamiltonian 
H = Hq + iJso° + Hso^ + if cm + -f^imp- As for the free electron part, Hq, the commutator 
reads 

Po.el^^W (Al) 

The spin-orbit contributions read 

ijk 

[Hso^, c] = e ^ e,,fc(V,i;(^))v4,afcC, (A2) 

ijk 

and the electromagnetic field contribution is 

[Hem. c] = -z- V AiV^C. (A3) 
m. ' 
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The commutator for the creation operator is given by [-f/^, c^] = — [i?, c]^ Equation (fT3l) in 
the operator form thus reads (s° = cV^c is the electron spin operator) 

dtS'' = TT^V ■ (cV« V c) - V e,,„(V,t;L")) V,(ctc) + - V A,\/,{c^a^c) 

(V.t;L^))(cV, v,c) + 2e5^ eijfcefcai(^i^io'')^i(c^'^«c). (A4) 

jj'fci ijkl 

(The contribution from the impurity scattering, i^imp, vanishes.) The first three terms on 
the right-hand side of Eq. f lA4p are written as divergence (choosing the electromagnetic 
vector potential as divergenceless, V ■ A = 0). From Eq. IA4t it is useful to define the spin 
current operators as 

(A5) 

where 

Taking the average, the equation of motion for the spin density operator, Eq. flA4p . results 
in the spin continuity equation, 

d,s- ^ -V ■ (jt^'" + + %: + X^'", (AT) 
where each terms are given by Eqs. (ITSi) (fTSl) . 



Appendix B: Calculation of spin relaxation torque 

In this section, we show details of calculation of spin relaxation torque. The first contri- 
bution in Fig. [Dreads 

kk' q ijkl mn^'y 

X Ai{k - k')i{k + k' + q)jk'Jtn (k' + ^^A'^i39k,Lo(r'ygk' ,u>gk'+q,Lo+n]'^ + cc. (Bl) 

Here guui represents the contour ordered electron Green's function with the wave vector k 
and the frequency oj {gku is a 2 x 2 diagonal matrix in the spin space), tr is the trace over 
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the spin, and [ ]^ represents the lesser components, and ( ); represents the averaging over 
the random spin-orbit impurities. Taking the lesser component, we obtain 

kk' q ijkl mnf5j 

X Ain{k - k'Uk + k' + q)jk'^kn {k' + I) ^ tT[(T pgla^gl,gl,+^] + c.c. (B2) 



The retarded and advanced electron Green's functions at zero frequency are represented by 
and g^, respectively. 

The leading contribution for small l/(eirr) is given by 

7T 777- ^ ^7T 

kk' Q ijkl nvnp'j 

X Aiflkik'jk'^knk'itT[a^gla^gl,gl,^g] + c.c. (B3) 
Expanding q in the Green's function, we obtain the leading contribution as 

kk' q ijklomnfi^ 



X 



Am,k'.k'^knkXqotT[apgla.,gl,{gl,Y] + c.c. (B4) 



Using the rotational symmetry, i.e., {kikj) = ^6ij and (kikjkikk) = ^{Sij5ki + 6ikSji + SuSjk) 
(the average here denotes the angular average), we obtain 

kk' q kll3-y 

X ^Aink\k')\-46k^qi + 6i^qk + Skiq^)tr[a^gla^gl,{gl,f] + c.c. (B5) 



We carry out the average over the impurity spin-orbit interaction as Eq. (jlOl) to obtain 
U = ^r.3oAso^^ I'^Y.Y.H ^^^.e-^'-'^e^'^'^iQ " A - q,A,/2)e{krtr[a,gla,gl,{gl,r] + c.c. 

kk' q 13-y 

(B6) 

The asymmetric part of the trace is calculated by use of {A{= ^{A^ + A_ + az{A^ — A_))) 
and B{= + B_ + a^iB^ — B_))) are any diagonal 2x2 matrices) 



tr 



[{a,3Aa^ - a^Aap)B] = -2iep^, ^ aA^B^,. (B7) 
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The result is 

fcfc' q a 

(B8) 

which is Eq. ([20]). 

Since the summation over k and k' is dominated by the imaginary part of the Green's 
functions, we use gl^^ = -^\gla? + o{^) and gi,^^M',~ay = ih^^l'^-J^ + 0(7^), and 



obtain the leading contribution as 

(B9) 



'so /It- "■so^so 2 

45 



Here spin-dependent density of states and the Fermi wavelength are represented by z/q- and 
kpo-, respectively. 

Another contribution to the relaxation torque, 7^o of Eq. (fT9|) . represented by the third 
diagram of Fig. [H reads 

fcp qQ ijk mnf}^ 

(BIO) 

As is seen, T^^'" is an odd function of k or k', and thus it vanishes. 
The second diagram in Fig {1] (we denote it as ST^o) reads 



^.■.^./3.e...e---e^-' {v^{l^ - k - q)vi2\k - k' + q)).^ 



kk' qQ ijk mn/3'y 

X Ai{q, n){k' -k- q)i{k' -k- q)n{k + A;')jtr[o-^5(fc,^cr^fi(fc/,<^+n]^- (Bll) 

Here, the wave vectors of the two Green's function, k and k' , are independent, and thus the 
magnitude of 57^" turns out to be smaller than T^l by the order of l/lepr). 

We therefore see that the dominant contribution to the relaxation torque arises from the 
first process in Fig. [H which is estimated to be Eq. (1B9[) (Eq. (!20|) ). 
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Appendix C: Calculation of the spin current 



1. Normal spin current 



Here we show first the derivation of the normal spin current without the vertex correction, 
i.e., the first term of Eq. (^^. The leading contribution of Eq. f lTSj) arises from the first 
diagram of Fig. El It reads 



■{o),z ■ e 



X 



kq j a 



+ f{ui — —]gt q n gt,q ,n — f ( u -\ — ] gl q n gl , q ,n 



(CI) 



Here the effect of spin-orbit interaction is included in the lifetime (selfenergy) of the Green's 
functions. Estimating j^^^'^ to linear order in and using f'{uj) = —6{ui) at low tempera- 
tures, we see obviously that the first contribution in the square bracket reduces to the first 
term of Eq. ( IMl) . The remaining contribution, which we call Sj^^^'^ , is calculated to linear 
order in Q as 



X 



kq - 3 

f (^) \9k~^,ui,a9k+i,uj,a ~ 9l-a.^^,a9k + %, 



~ '2^' l9k-^,LU,a9k+§,LU,a + 9k-§,uj,a9k+§,ui,o] 

-{9k-!i,u>,a)^9l+!l^uj,a + S'fc-f (S'fc+f 



(C2) 



On the other hand, the contribution from the gauge field, j^- in Eq. (1151) . reads 

, t) = iSa,z—Ai{x, t) f^y] ag^^^ 
m I In ^ — ^ 

^ ka 

m J 2,71 ^ — ^ 



■Am/ ^ 
JsJ (^v 



(C3) 



This contribution is proportional to A, and so is of the order of E/Q. One can easily see 
that this contribution cancels with the first term of Eq. flC2p estimated at g = 0. We 
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therefore obtain 



Js,t ' Js,t 



X 



kq j a 



n 



+ i^) [(fi'fc-f 9k+^,uj,a ~ 9k^^,uj,ai9k+^,uj,a) 

-(^fc-f ,a.,a)'^fe+f + 9l-^,^j9l+^,^J^] + 0{q) 



(C4) 



Since we consider a slowly varying gauge field, it is sufficient for our purpose to estimate 
this expression at g = in the Green's functions (the adiabatic limit). The result is 



Js.i ' Js,t 



kq j a 



x{-^n^)]M,.,.r + i9u.)% 



(C5) 



which turns out after fc-summation to be negligibly small (smaller by a factor of l/{eFT) 
than the dominant contribution). We therefore obtain Eq. (1241) as the dominant normal 
spin current at small Q and q. 

2. Calculation of the third diagram in Fig. [2] 



The contribution from the third diagram in Fig. [21 which arises from the normal current 
and the spin-orbit interaction including the gauge field, reads 



J sA 



kk' q 



~ ~ I) j^^^^^"9k-lc.'T,39k'^lco(^-,9k+^ 



(C6) 



p=k'—k 

Taking the lesser component, the dominant contribution at small Q and q turns out to be 



.(3),a rdu r dQ 

" m / 2^ / 2^ 



kk' q 



(C7) 
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By use of the averaging over the spin-orbit potential diagonal in the spin space, Eq. (fTOjl . 
we obtain 



.(3),a _ . • 



•J UUI n 



kk' q 

X ti[(j''glap{gl. + gl,)(y^gl]p^^. (C8) 

We therefore see by noting that 'Yl,k'i9k' +fi'fc') — Oi^u I {^f't)) that the contribution is 
negligibly small compared with the main contribution, Eq. 



3. Anomalous spin current 

The anomalous spin current arising from the spin-orbit interaction, defined in Eq. ( |T6|) 
and shown in Fig. [2] (the fourth and fifth diagrams), is calculated in a similar manner as 
that for the relaxation torque. The dominant contribution arises from the fourth diagram 
of Fig. [21 It reads 



•so, a -(4), a 



kk' q 

X Ai{k - k')jk'^kn (k' + ^'^i[gk,oj(r-~igk' ,u^9k'+q,ui+n\^ ■ (C9) 

Taking the lesser component and the average over the spin-orbit potential, and by noting 
that only 7 = 2; contribution is finite because of the trace over spin, we obtain 

= ^^soAso'- E E e.,.e_,e-'^-e^™A,fi(A: - k'),k'^k^ {k' + |) tT[gla,gl,gl,^^] 

^ 1,1.1 n ' 

(CIO) 



kk' q 



The leading contribution at small q is then obtained as 

-i ,96 

kk 



= ir:nsoKo^—Mx^i)^^ijz{kjf{K)hY[gla^gl,gl] 



-f nsoAso^^ii(a;, t) (T{kFaY{VafT^ (« = X, y) 



(Cll) 



U = ^ 



where we used the rotational symmetry in k space. This expression is Eq. ( 135|) . 
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The last diagram of Fig. |2] reads 

Saz-^nsoXso^etjzekizA ^(/c - k')j{k - h')}M\g\(j ^.gW 



.(5),a 
? • 



(C12) 



fcfe' 



This contribution turns out to be smaller than Eq. (IClip by the order of (cft) ^• 



Appendix D: Calculation of asymmetric /3 torque 

The torque generated by the spin-orbit interaction and uniform electric field when the 
spin structure is inhomogeneous was calculated in Ref. 10|. There, the angular averaging 
over the electron wave vectors was carried out assuming rotational symmetry on averaging 
the spin-orbit interaction, resulting in a symmetric torque called the /3 term. Here we carry 
out the calculation without such approximation, and show that asymmetric /3 term then 
arises. The torque reads (see Eqs. (19)(21)(35) of Ref. |lO| ) 



7rm 



2 '^so'^so 



Iiul3-yij kk' 

X tT[(apgl,a^ + (T^gl,crp){\gl\^asgl + gl(T5\gl\% 



(Dl) 



where R^^, = 2m^m,^ — 6^^, denotes a rotation matrix for the gauge transformation and 
Aj = (m X VjinY is the spin gauge field (m = (sin | cos 0, sin | sin 0, cos |), where 9 and 
are polar coordinates representing the magnetization direction) 10[. In Ref. 10|, the term 
[{k X arising from the spin-orbit interaction was averaged independent of the other 

factor, kikj, and thus the torque was proportional to (kikj) = ^6ij (Eq. (37) of Ref. 10|). 
Here we carry out the averaging without assuming special symmetry. The averaging over 
the wave vectors in Eq. (]D1|) then reads 



{k,k,[{k X k%]^) = -k\k')\26,,-6,,6,,). 



(D2) 
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The last term gives rise to asymmetric torques. Evaluating the trace in the spin space, the 
torque reads 

8e 



T" = n 



^uP^ ^ ^ kk' 



cr' 



(D3) 



Symmetric contribution (contribution from the E ■ term) to the first line in the square 



bracket was calculated in Ref. 



10| by use of 



W-^M/3^-7^i - ^Piz^&z) = 2(v4f + n'^A^) = (n X V^n)", (D4) 

where n denotes the magnetization direction. ( The second and third terms in the square 
bracket in Eq. fIDSp lead to a correction to the spin Berry phase (proportional to A^) and the 
standard spin transfer torque (proportional to Vn"), respectively. ) As for the asymmetric 
contribution (terms proportional to ^Ei,Al in Eq. (IDSP ). summation over u cannot taken 
independent of the electric field, resulting in complicated torques. We here look into only 
the correction to the /3 term (the first line in the square bracket), which is obtained as 
(neglecting irrelevant terms) 



(D5) 



The result of the /3 term is thus 



7;(/3)." ~ _[/3n X {j ■ V)n + 5(3n x (j||V||)n] 



(D6) 



where j\\ = n ■ j, Vy = n ■ V, 

8e 



/3 



(D7) 



kk' aa' 

{as is the Boltzmann conductivity) and 6P//3 = —1/5 in the present model. Thus the /3 
term arising from the spin-orbit interaction has asymmetric contribution, in contrast to the 
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one arising from the random magnetic impurities 
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